Critical Phenomena in a Highly Constrained Classical Spin System: 
Neel Ordering from the Coulomb Phase 
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Many classical, geometrically frustrated antiferromagnets have macroscopically degenerate ground 
states. In a class of three-dimensional systems, the set of degenerate ground states has power-law 
correlations and is an example of a Coulomb phase. We investigate Neel ordering from such a 
Coulomb phase, induced by weak additional interactions that lift the degeneracy. We show that the 
critical point belongs to a universality class that is different from the one for the equivalent transition 
out of the paramagnetic phase, and that it is characterised by effective long-range interactions; 
alternatively, ordering may be discontinuous. We suggest that a transition of this type may be 
realised by applying uniaxial stress to a pyrochlore antiferromagnet. 

PACS numbers: 64.60.Cn 75.10.Hk 05.70.Jk 



The statistical physics of systems that combine local 
constraints with macroscopic degeneracy has been of en- 
during interest in several different contexts. Examples 
include models for geometrically frustrated antiferromag- 
nets [l[ and for proton disorder in ice J3] , and dimer cover- 
ing problems on a variety of lattices [3| ■ Local constraints 
arise from energy minimisation when ground states of a 
Hamiltonian are considered, or are rules that must be sat- 
isfied by all allowed configurations, as for closed-packed 
dimers on a lattice. Macroscopic degeneracy appears if 
there are extensively many ways to satisfy these con- 
straints. Despite the degeneracy, correlation functions, 
averaged over states of this kind, need not be short-range. 
Indeed, it has recently been appreciated, both for dimer 
coverings of three-dimensional bipartite lattices, and for 
ground states of certain geometrically frustrated mag- 
nets, that pair correlations at large separation have a 
characteristic, algebraic dipolar form 0,(1,0. This be- 
haviour is captured by the Coulomb phase of a theory 
in which the local constraints of lattice models become 
the requirement that a continuum field is solcnoidal. The 
Coulomb phase is stable against all small perturbations 
that are allowed by symmetry, but a sufficiently large 
perturbation may induce long range order. 

In this paper we consider a Neel ordering transition 
that takes place from the Coulomb phase of a geomet- 
rically frustrated antiferromagnet, asking how it differs 
from the conventional Neel transition out of the param- 
agnetic phase. Our interest in such a transition stems in 
part from the fact that, by design, it evades a naive ap- 
plication of Landau-Ginzburg- Wilson theory, since a de- 
scription based in the standard way on symmetries and a 
local order parameter for the Neel phase will not repro- 
duce the correlations of the Coulomb phase 0] • Similar 
considerations were also a motivation for an intriguing re- 
cent study of another ordering transition from a Coulomb 
)hase, involving close-packed dimers on a cubic lattice 
|| • In that case, high-precision Monte Carlo simulations 
indicate critical behaviour different from that expected at 



a conventional critical point with the same symmetries, 
but an analytical understanding of the transition has yet 
to be developed. By contrast, we show in the following 
how the Neel ordering transition from the Coulomb phase 
may be brought within the standard theoretical frame- 
work used to describe critical phenomena. More broadly, 
in light of the intense current interest in quantum phase 
transitions that lie beyond the scope of a conventional 
Landau treatment f9|, it is natural to examine ordering 
from the Coulomb phase as a classical counterpart. 





FIG. 1: Neel order induced by uniaxial compression in the 
building blocks of geometrically frustrated antiferromagnets. 
Left: a tetrahedron from the pyrochlore lattice strained along 
[100]. Right: an octahedron from a lattice of corner-sharing 
octahedra strained along [111]. Full lines denote antiferro- 
magnetic exchange interactions of strength J, and dashed 
lines indicate ones of strength J — A. Arrows represent spins. 

The class of models we are concerned with is as 
follows. We study classical antiferromagnets on geo- 
metrically frustrated, three-dimensional lattices that are 
known to support a Coulomb phase in the zero temper- 
ature limit when all non-zero exchange interactions are 
the same. These systems meet three requirements: they 
have macroscopically degenerate ground states, do not 
display order by disorder, and have a bisimplex lattice. 
Such lattices consist of corner-sharing frustrated units, or 
simplices, which are arranged in a bipartite fashion, with 
exchange interactions linking all pairs of spins within 
each simplex. We focus on two examples: the Heisen- 
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berg model on the pyrochlore lattice, built from corner- 
sharing tetrahedra, which is approximately realised in a 
number of materials; and a system of rt-component spins 
on a lattice of corner-sharing octahedra, which is a con- 
venient theoretical construct. The local constraints sat- 
isfied in ground states of these models are that the total 
spin of each simplex is zero [Io| . We induce order using 
a change in the relative strength of exchange interactions 
which can be produced by uniaxial compression via mag- 
netoelastic coupling, as illustrated in Fig. [TJ Exchange 
energies Jij within each simplex then take the two values 
J and J — A, and our focus is the regime J ^> A > 0. 
The Hamiltonian for both models has the form 



(1) 



where Si is an n-component unit vector. In each case, 
positive A removes the macroscopic ground state de- 
generacy of the unstrained system, selecting from these 
ground states the two-sublattice Neel states of Fig. [TJ 
in which neighbouring pairs of spins coupled by the 
stronger interaction J are antiferromagnetically aligned, 
while frustration forces those pairs coupled by the weaker 
interaction J — A to have ferromagnetic alignment. Con- 
sidering the effect of thermal fluctuations, one expects 
Neel order below a critical temperature T c ~ A, and a 
disordered phase at higher temperatures. For temper- 
atures T <C J, thermal fluctuations are predominantly 
within the ground state manifold of the unstrained sys- 
tem, and in the limit J/A — » oo the transition at T c is 
between a Neel phase and a Coulomb phase. 

Several examples are known of magnetic ordering asso- 
ciated with lattice strain in frustrated magnets, although 
we are not aware of a direct realisation of the transition 
we discuss. In particular, Neel order has been induced 
by uniaxial stress in the pyrochlore Tb2Ti207 while 
it is accompanied by spontaneous strain in the spinel 
ZnCr204 [12J. Unfortunately, Tb2Ti2C>7 is a material 
that is at present poorly understood, and the transition 
in ZnCi'204 is strongly first order. We hope that this 
paper will stimulate further experiments in the area. 

We set out below several complementary approaches to 
understanding ordering in constrained systems. Starting 
from the continuum description of the Coulomb phase, 
Neel order amounts to condensation of the solenoidal 
field. When this transition is continuous, we show that it 
has the same critical behaviour as for an unconstrained 
system in which there are uniaxial dipolar interactions. 
From classic studies [H, HI 15| . it is known that long- 
range interactions of this type lead to an upper critical 
dimension of three, and we confirm this for the lattice 
models in the large-n limit. We supplement this analysis 
with Monte Carlo simulations. We find that the transi- 
tion in a Heisenberg model is first order, but show that 
a critical point can be reached by applying a staggered 
field. 



We start by recalling the mapping of Refs. [B| and @, 
from ground states at A = to solenoidal fields. The 
mapping uses the underlying simplex lattice for the spin 
system, which is a diamond lattice for the pyrochlore an- 
tiferromagnet, and a simple cubic lattice for the octahe- 
dral antiferromagnet. In each case, spins Sj are located 
at mid-points of bonds i on the simplex lattice. Unit vec- 
tors §i are defined on these bonds, with the orientation 
convention that they are all directed from a given simplex 
sublattice towards the other one. A spin configuration is 
represented using n flavours of a three-component field, 
defined from the spin components Sf by B™ = Sf&i. 
The ground state constraint, that the total spin of each 
simplex is zero, translates into the condition that these 
fields have lattice divergence zero. Introducing contin- 
uum coordinates r and coarse-grained fields B a (r) satis- 
fying V • B Q (r) = 0, the entropic weight on ground states 
is postulated to be proportional to e~ H °, where 



d 3 r^|B Q (r)| 2 



(2) 



The stiffness k is the single parameter of the theory: it 
simply controls the amplitude of correlations. 

Under the same mapping, the Neel states of Fig. [1] are 
represented by a uniform, non-zero B Q (r). The orien- 
tation of this field in real space is along [100] for the 
pyrochlore antiferromagnet, and along [111] for the octa- 
hedral antiferromagnet, while the orientation in flavour 
space reflects the O(n) symmetry-breaking of the Neel 
order. 

In order to describe a phase transition between the 
Coulomb and Neel phases, we incorporate the effect of 
non-zero A into the continuum theory of Eq. ([2]) via a 
temperature-dependent control parameter i, with t = 1 
for T > A and * < for T < A. Separating B Q (r) into 
components -Bfj*(r) and B" (r), respectively parallel and 
perpendicular to the preferred direction, we replace Ho 
by H2 + Ha , where 

H 2 = ^ / d 3 r^{|B^(r)| 2 +i|5f(r)| 2 + |VxB«(r)| 2 } 

H 4 = u / d 3 r{^|B-(r)| 2 } 2 . (3) 

a 

At t = 1, Eq. ([3|) is simply Ho supplemented by gradient 
and quartic terms that are irrelevant in the scaling sense 
in the Coulomb phase. On reducing t a transition occurs 
(at t = within mean field theory) to an ordered phase in 
which (B^(r)) ^ 0. The quartic term H4 is required for 
stability in the ordered phase, while the gradient term is 
necessary to suppress short-wavelength fluctuations close 
to the critical point. 

The solenoidal constraint can be re-expressed as a long- 
range interaction. To this end, write the fields B Q (r) 
in terms of vector potentials A a (r) by setting B a (r) = 
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V x A Q (r ) , and choose the Coulomb gauge V • A a (r ) = . 
Then TL2 is diagonalised by Fourier transform and by 
rotation in the space of the three components of A". 
In a fixed gauge only two of these components fluctuate 
independently: taking their amplitudes at wavevector q 
to be a"(q) and ci^q), the diagonalised form is 

U * = \ I [H +9 4 K(-q)o?(q) 

+ [<Z 2 + g 4 K(-qK(q)}. (4) 

To emphasise the equivalence to a system with uniaxial 
dipolar interactions, we rescale the fields by a factor q, 
writing <pf(q) — qaf(q) for I — 1,2; and since <£>f (q) is 
non-critical, we omit it, obtaining finally 



n-> 



^E 



(l-t)4+i + g 2 



^r(-qK(q)- 

(5) 



In this form, Ti. 2 has been studied extensively as a 
model for uniaxial ferroelectrics and ferromagnets, with 
the term Q 2 /q , non- analytic at q = 0, arising from dipo- 
lar interactions. Ordering in these systems is celebrated 
as an example of a phase transition with an upper critical 
dimension of three, at which critical behaviour is almost 
mean field, but with logarithmic corrections 13, 14, la ]. 
We conclude that the same behaviour will appear at the 
Neel ordering transitions of the constrained spin systems 
discussed here, a demonstration of the controlling influ- 
ence that the correlations of the Coulomb phase have on 
critical phenomena. 

We next turn to lattice models, and first present an 
alternative way of viewing the equivalence between the 
constrained spin systems we study and ones with uni- 
axial dipolar interactions. Consider the the spectrum of 
the interaction matrix Jy appearing in Eq. ([IJ. It is 
block-diagonalised by Fourier transform with wavevector 
k. In the class of models we consider, at A = its lowest 
branch is dispersionless, and for bisimplex lattices this 
branch is degenerate with dispersive branches at points 
ko in the Brillouin zone. A natural approach to statisti- 
cal mechanics at T <C J is via projection onto these flat 
bands. The projection operator P is non-analytic in k at 
kn and as a result is power-law in real space at large dis- 
tances 16]. For three-dimensional, bisimplex frustrated 
antifcrromagnets, its long-distance form is dipolar 17], 

The dipolar nature of P, applied to certain pyrochlore 
magnets (spin-ice materials |l8|) turns out to have the 
remarkable consequence that dipolar interactions are es- 
sentially equivalent to nearest-neighbour ones, within the 
set of low-energy states allowed by local constraints, a 
feature termed projective equivalence '. 
an opposite path: we show that projection of Jy 
A > turns nearest-neighbour interactions into long- 
range ones. 



Here we take 
for 



On favouring Neel order by setting A > 0, degen- 
eracy in the spectrum of Jy is lifted, leaving isolated 
minima. In this language, the models of Fig. Q] are en- 
gineered so that these minima occur at the degeneracy 
points ko- Then for A <g. J the interaction matrix in- 
herits the singularities of the projection operator. To be 
definite, we now focus on the octahedral lattice, for which 
ko = (7r, 7r, 7r); we have also obtained equivalent results 
for the pyrochlore lattice. Introducing q = k k?,-, the 
projected interaction has a lowest eigenvalue that varies 
with q as 

£(q) = £ + 3A4 + 0{q 2 ) , (6) 



where q\\ is the component of q in the [111] direction and 
£0 is a constant. Spin configurations that are thermally 
allowed for T / J — > are invariant under projection: the 
presence in Eq. © of the singular term q 2 / q 2 is therefore 
a key result, indicating effective uniaxial dipolar interac- 
tions within these configurations. 

With T/J small but non-zero, the power-law correla- 
tions of the Coulomb phase are cut off at a large distance 
£c- I n these circumstances, uniaxial dipolar interactions 
will control critical phenomena provided the correlation 
length £ is smaller than £c , but with a crossover to con- 
ventional behaviour sufficiently close to T c . 

To study the phase transition in the lattice model, we 
use the large-n limit, extending the treatment of geomet- 
rically frustrated systems described in [B|, , and recov- 
ering the critical behaviour known for uniaxial dipolar 
systems at n — > 00 13j, [lJ, |l5| . As a function of reduced 



temperature i r , we find that the staggered susceptibility 
X diverges as x ~ |* r | 1 • ln(l/|t r |), and the sublattice 
magnetisation vanishes asm~ 1 | 1 / 2 . 

To examine these ideas further, we use Monte Carlo 
simulations. We study the pyrochlore Heisenberg anti- 
ferromagnet, treated extensively at A — in previous 
work [lOj. Crucially, the absence of energy barriers in 
this model at A = makes it possible to equilibrate 
simulations at temperatures as low as T = 1CP 4 J, per- 
mitting access to the regime we are concerned with. We 
focus on the two choices: A = 0.1J and A = 10~ 3 J, us- 
ing systems of N = 4L 3 spins, with N < 4000 (L < 10) 
and employing parallel tempering [2(]| . 

For A = 0.1 J, we find (data not shown) a continuous 
transition at T c = 0.267J and critical behaviour con- 
sistent with the standard, short-range universality class, 
which is unsurprising since £c is not large at this tem- 
perature (£ c oc (J/T) 1 / 2 for n > 2 and T < J). For 
A = 10~ 3 J the data illustrated in Fig. [2] indicate a 
strongly first-order transition. The Neel order parame- 
ter shows a rapid change with temperature, and this be- 
comes more step-like with increasing system size. More- 
over, close to the transition the probability distribution 
for the energy develops two peaks, a feature characteris- 
tic of a first-order transition [2l| . 
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FIG. 2: Evidence for a first-order transition from simulations 
with A = 1CT 3 J. Main panel: m vs. T/J for L = 6 (□), 
8 (■), and 10 (o). Inset: probability distribution p(E/J) 
for the energy E at the estimated transition temperature, 
T/J = 4.6 x 10~ 3 , for L = 6. 



There is an obvious interest in finding a version of this 
ordering transition which is continuous even at A < J, 
as a realisation of the uniaxial critical behaviour we have 
discussed. It is possible that the order of the transition 
could be altered by a change in microscopic details of 
the system. Here, however, we follow an alternative and 
more automatic route to a critical point, by applying a 
staggered field: we add to Eq. |T]) the term —hJ2i n iSi> 
with rii = ±1 taking opposite signs on the two Neel sub- 
lattices. For h -C A, a first-order transition survives, 
between a high-temperature phase with small sublattice 
magnetisation m and a low-temperature one with large 
m. Increasing h, the discontinuity in m at the transition 
is reduced, and above a critical field h c ~ A there is no 
thermally driven phase transition. The critical end point 
at h = h c is the continuous transition we seek. Evidence 
for its existence is presented in Fig[3J A precise numer- 
ical study of critical behaviour at h — h c for A/J « 1 
would be challenging and we do not attempt it here, but 
we believe the theoretical case for uniaxial dipolar critical 
behaviour at this transition is compelling. We note that, 
while the staggered field used in our simulations cannot 
be realised in a Heisenberg system, it would arise from a 
uniform external field in spin- ice materials (l8j . 

In summary, we have shown that ordering in highly 
constrained spin systems can be quite different from its 
counterpart in conventional systems. We suggest that 
experiments on ordering in uniaxially strained geometri- 
cally frustrated magnets should be of great interest. 
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sions. This work was supported by EPSRC Grant No. 
GR/R83712/01. Part of it was done while JTC was a 
visitor at KITP Santa Barbara, supported by NSF Grant 
No. PHY99-07949. 
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FIG. 3: Transition in a staggered field with A = 10 -3 J. Left 
inset: Schematic phase diagram, showing first-order transi- 
tion line terminating at a critical end point. Main panel: 
smooth variation of m with T at h — 2 x 10 _3 J, which 
is greater than h c , for L = 6 (o) and L = 8 (•). Right 
inset: probability distribution p(E/J) for the energy E at 
T/J = 8.0 x 10~ 3 , 8.3 x 10" 3 and 8.6 x 10~ 3 , showing single 
peaks. 
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